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\Q ■ ABSTRACT 

Aims. We present the basic equations and sample solutions for the steady-state global transonic 
solutions of the advection-dominated accretion flows (ADAFs) with causal viscosity prescrip- 
tion. The procedures of the stable numerical calculations and all explicit formula to obtain the 
2 . solutions of ADAFs are also presented. 

Methods. We solve the transonic solutions of ADAFs by using the explicit numerical integra- 

C3 . 

tions, such as the Runge-Kutta method. In this calculation method, we firstly solve the physical 
values at the sonic radius where L'Hopital's rule is used. Then, we numerically solve the coupled 
^ ' differential equations of the radial velocity, the angular momentum and the sound speed from the 

wn 

■ Results. By the calculation procedures presented in this paper, we can cover the all parameter 



spaces of the transonic solutions of ADAFs. Sample transonic solutions for ADAF-thick disk and 
ADAF-thin disk are presented. The explicit formula for the analytical expansion around the sin- 



o . 

, gular points, the sonic point and the viscous point, are presented. If we set the diffusion timescale 



to be null, the formalism in this study become the formalism of the acausal viscosity which is 
usually used in the past study for the calculation of the ADAF structure. 

Key words, accretion: accretion disks — black hole physics — hydrodynamics 

1. Introduction 

Advection dominated accretion flows (ADAFs) are one of the fundamental solutions of the accre- 
tion flows into the black holes in the low-luminous systems, such as, low-luminosity active galactic 
nuclei. In past studies, the several methods are used to obtain the transonic solutions of ADAFs, 
e.g., the relaxation scheme (e.g. Narayan, Kato & Honma 1997), the implicit integration scheme 
(e.g. Kato, Fukue & Mineshige 1998), the outward integration scheme with the inner boundary 
condition (Becker & Le 2003) or the inward or outward integration scheme with the special con- 
dition at and inside the sonic radius (Lu, Gu & Yuan 1999). All these works assumes the diffusive 
viscosity with infinite diffusive velocity. 

In this paper, we give the method to solve the global transonic solutions of ADAFs with 
causal viscosity prescriptions which assume the finite diffusive velocity proposed by Papaloizou 
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& Szuszkiewicz (1994), which can become the alternative method to solve the global solutions. 
Gammie & Popham (1998) calculate the ADAF structure in the Kerr metric with the causal vis- 
cousity prescription by using the relaxation scheme. In the calculations with the causal viscosity 
prescription, the boundary condition at the viscous point is used, and the boundary condition at the 
horizon is not required (see, also §5). If we set the diffusion timescale to be null, the basic equations 
in this study become the basic equations for the case of the acausal viscosity which is usually used 
in the past study. 

After we present the basic equations for the ADAF with the causal viscosity prescription in §2, 
the boundary conditions for the global transonic solutions of ADAFs with the causal viscosity are 
given in §3. The numerical procedures of the method to solve the transonic solutions of ADAFs 
and the sample solutions are presented in §4. Discussion and conclusions are given in §5 and §6, 
respectively. 



2. Basic Equations of ADAF with Causal Viscosity 

The basic equations used in this paper is same as those in Narayan, Kato & Honma (1997) except 
that we use the causal viscosity prescription. The mass conservation is written as, 

- AnrHpv r = M, (1) 

where H is the half thickness of the accretion disk, p is the gas density, v r is the radial velocity 
and M is the mass accretion rate. Here, we assume v r < and H - (5/2) 1 ^ 2 a J /Q^ where a s is the 
sound speed defined as a s = (plp) l l 2 and p is the pressure. 
The radial momentum conservation is described as, 

dVr <cV n2^ 1 d P ™ 
v r — =(Q -Q K )r- - —, (2) 

dr p dr 

where £2 is the angular velocity of the accretion flow and is the Keplerian angular velocity 
calculated as Q^- = l/[r 1/2 (r - 2)]. Here, we set GM = 1 and use the Paczynskii-Wiita potential. 
The energy equation is written as, 

pv r da] 2 dp dQ 

—l^r-- a s Vr dr= PVrie - j) ^> (3) 
where y is the ratio of specific heats of the accreting gas, €{= Qr 2 ) is the angular momentum, and 
j is the integration constant meaning the specific angular momentum at some radius. 

For the case of the diffusive viscosity with the infinite diffusion velocity, the shear stress t r $ is 
written as 

t r<p = -vpr*-^. (4) 
Then, the angular momentum equation for the case of the acausal viscosity is written as, 

dQ. _ Vrd-j) 

J 2 ' ^ ' 

dr vr^ 

where v is the kinematic coefficient of viscosity, and, here, j is the specific angular momentum per 
unit mass where dQ./dr = 0, if such point exists. By assuming the a-viscosity we set v = aa^/O.^. 
From Eq. (O, the equation for the angular momentum is written as 
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In the same way as Papaloizou & Szuszkiewicz (1994) and Gammie & Popham (1998), we 
include the effects of the diffusive viscosity with the finite diffusion velocity. By using the causal 
viscosity, we can derive the angular momentum equation with the diffusive viscosity with the dif- 
fusion velocity a v written as, 



<M__Ny_ 

dr ~ D Y ' 
where 



(7) 



<n 1 2,2 Kr 21 v r (€ - j) ( 2v r T v 
D v = l-v r /a y , N v = — + 1 , (8) 



Here, r v is the relaxation timescale. For simplicity, we set t„ and a v as r v = Q^. 1 and a v = (v/Vy) 1 ' 2 , 
respectively. Now, the shear stress for the causal viscosity is written as 



_ -yp 



[ v 2 \ 2 dQ. 2£v 2 
dr r a 2 



(9) 



1 - 2v r T v /r 

When t v - 0, i.e., a v - oo, Eq. (0 and (O for the causal viscosity become the equations for the 
acausal viscosity, Eq. © and Even in the point where d£l/dr = 0, the shear stress t r< f, in general 
do not become null due to effects of the causal viscosity. Since, the effects of the causal viscosity 
are coupled with the radial velocity as v 2 /a 2 and v,r v in Eqs. (0 and ©, the effects of the causal 
viscosity is effective in the inner region of the accretion flow where the absolute value of the radial 
velocity is large. 



3. Boundary Conditions and Physical Values at the Sonic Radius 



Since the accretion flows plunge into the black hole supersonically after a sonic transition, the 
equations are singular at the sonic radius. In order to see this, from Eqs. (Q]), and (0, we can 
derive the equation for dv r /dr as, 



dv r 
~dr~ 

where 



N 



(10) 



D = v r - 
Q 2 - 



2y \a] 
y + 1 ) v r 



N 



2y \ 2 dhW*/r) 



y+l 



dr 



(y-l\ (e-j) (d( 2£\ 
\y+ 1 ) r 2 \dr r )' 



Here, Eq. (O is inserted to dljdr and d\n{Q.K I r) I dr 
the equation for da jdr is written as, 

fy-V 



(11) 



-(5r - 6)/[2r(r - 2)]. On the other hand, 



da s 
dr 



\y + 1 



a 

Vr 



s dv r dln(a K /r) (t - j) Idt 2£\ 
r dr s dr a s r 2 \dr r j 



(12) 



In order to pass the sonic radius, r s , smoothly, two boundary conditions at r — r s are given as 
T> — N — at r — r s . For given a s>s (> 0) which is the sound speed at the sonic radius, the radial 
velocity v r , s at r — r s is calculated from T) = as v r , s = —o s , s [2y "/ (y + l)] 1/2 - From the condition 



4 R. Takahashi: Advection-Dominated Accretion Flows with Causal Viscosity 

N — 0, the angular momentum i s at r — r s is calculated as, l s - [be + (b 2 - aece) l ^ 2 ]/ae where 



at = 1 - 
be = 



7-1 
y + 1 

7 / 1 - 7 



rv r , 
2+ — 11- 

V 



2v r T r 



Cf = - 



2 \1 +y 
-2 + ± 

y - 1 \ rv,.;' 2 





I 1 - 2 ";") 




V 





1 - 



2v r T, 



y + 1 / vD v 

2y 2 ^]n(Qjc/r) 
-arr 



-Cl 2 K r* 



(13) 



(14) 



(15) 



y + 1 dr 

Chakrabarti introduce very clever method to calculate the value of dv r jdr at r — r s (e.g. 
Chakrabarti 1990, 1996). We extend their method to calculate the transonic solutions of the ADAFs 
with causal viscosity. The value of dv r /dr at r — r s is calculated by applying L'Hopital's rule as 
follows. To apply L'Hopital's rule to Eq. (f2j) at r — r s , we introduce new variables x = r - r s , 
y = v r - v r>! , z — t -i s and w - a s — a 3iS . The functions D and N are expanded near the sonic point 
as, 



D : 

N : 
where 

\dx 

m 



idN_ 

\dy 
(dN 
\ dw 



<m 

dx 
dN 
dx 



x + 



X + 



3D 
~dz 
1 



3D 
~dy~ 

9y) s 



y + 
I y + 



3D 
~~dz 
dN 
~dz 



z + 



z + 



= 0, 



dD\ 



= 2, 



3D 
dw 
dN 
dw 



3D 
dw 



w, 



= 2. 



(16) 
(17) 



2y 
'y+ r 



= 



(r-2) 3 



3f 2 



2y \ 2 d 2 ln(£l K /r) 

a. 



y + 1 



dr 2 



y-l\6€(€-j) 



y + 1 
3N 
~~dz 



r 4 

2t 
z 3 



y-l\2(2t-j) 



-4y\ d\n(Q K /r) 



y + 1 



dr 



(18) 



Here, d 2 \n(Cl K /r)/dr 2 = (5r 2 - 12r+ 12)/[2r 2 (r-2) 2 ], andz and w are described as z = (d£/dr) s x 
and w = (da s /dr) s x + (da s /dv r ) s y where 



( da s 
\~dr~ 

( da s 



7 - 1 \ a s 



y + 1 / Vr 

y-i 

y + 1 



a dln(n K /r) | (i-j) 
s dr a,r 2 



21 
r 



Finally, we obtain the quadratic equations of dv r /dr at r — r s as, (dv r /dr) s = a + {a 2 - b) 1 ^ 2 which 
are the solutions of the equation (dv r /dr) 2 - 2a(dv r /dr) s + b — where 

\dyj s \dw) s \dv r )s \Sx) s \dz) s \dr) s \dw) s \dr) s 



a = 



b = 



2 [(fl?) s + (fl*)s(^v)J 

\ Ox Is \ dz Js\ oris V ow Is \ Sr ) s 

(¥) +{'¥) (¥) 

\dyJs \9wJ s \dv r J s 



(19) 
(20) 
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In addition to the boundary conditions T> — N — 0, when the causal viscosity prescription is 
adopted, the accretion flows smoothly pass the viscous radius, r v , where the condition T> v = is 
satisfied. In this case, the two boundary conditions at the viscous radius, r v , are given as D v — N — 
at r = r v . The method presented in this study uses the boundary conditions at the sonic radius, r s , 
and the viscous radius, r v , and do not use the boundary conditions in the outer region r > r y and 
the inner region r/, < r < r s , where r/, is the radius for the event horizon. 



4. Calculation Method and Sample Solutions 

Based on the formula presented in the previous sections, we can calculate the global transonic 
solutions for ADAFs with the causal viscosity. The calculation procedures to obtain the global 
transonic solutions for ADAFs are as follows: 

1. First, we tentatively choose some value of a ss for given values of r s and j, and calculate v ns , 
( s , (dv r /dr) s , (d£/dr) s and (da s /dr) s by using Eqs. (0, CHI and <fl2t , 

2. Next, we solve the solutions in the range r s < r < r v . In order to do this, we solve Eqs. (0, (TT0T > 
and ( fT2T > from the sonic point to the viscous point by using, e.g., the Runge-Kutta algorithm. 
Usually, for the initially selected value of a s>s , the calculated solution does not pass the viscous 
point where two boundary conditions D v — N v — are satisfied. In such case, we return to step 
1 and again choose the different values of a SyS for given values of r s and j. After repeating these 
procedures, we can determine the value a M which gives the solution satisfying the boundary 
conditions D = N = at r = r s and D v = N v = at r — r v . 

3. After solving the solutions in r s < r < r v , we solve Eqs. (0, ( TT0T > and ( TTZI ) in the range r v < r by 
using the values of a 3iS for given values of r s and j by using, e.g., the Runge-Kutta algorithm. 

4. Finally, we solve Eqs. (0, (fTOb and ( fT2T > in the range r/, < r < r s by using the values of a s>s for 
given values of r s and j by using, e.g., the Runge-Kutta algorithm. 

The third step and the fourth step can be interchanged. By this procedure, the transonic solutions 
are obtained for given values of r s and j. At the third step, the physical values, such as d£/dr, 
at the viscous radius r = r v are used if necessary to obtain the solutions. Usually, the numerical 
calculation of the coupled equations are naturally performed through the viscous point, and the cal- 
culations in the third steps are successively done after passing through the viscous radius. However, 
in some cases, the numerical calculation is stopped at the viscous radius. In such case, we evaluate 
the physical values at the viscous point by using the L'Hopital's rule which give the value dE/dr at 
the viscous radius r = r Y , and then we calculate the coupled equations for the global solution from 
the viscous radius to the outward direction r > r v . In the same way as the calculation of (dv r /dr) s 
as described in the last section, the value (d£/dr) r=r< . is calculated as (d(/dr) r=rv - a v ± (eft - b v ) 111 
where 

1 {ldN v \ (dN v \ (da s \ (da A (dD v 



2 \\ dz )„ \ dw M d( / \ dr / \ dx 



dD Y \ ( dD v \ l da s 
dy )y \ dw J v \dv r 

dDy \ tdD V \ IdOs 

dz )y + \dw)\d{ 
dN v \ t dN v \ (dc^ 
dx /., \ dw /„ I dr 



(21) 
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dy 

dDy 

~dT 



dN v 
dw 

dw 



da s 
dv r 

da s 
~~dt 



K 

D v 



(22) 



Here, 

\ dx 
( OK 
[ dx 
dN v 
Oy 
dN v 
dz 
dN v 
dw 
da s 
dr 
da s 
dv r 



dD 
~dz 

_K_ 
r 

(i-j) 



= 0, 



dD, 
~dy~ 

({-]) d(£l K r) 

rv r dr 
v r (Q/j-r - 2v r ) 



a x l 2 a 



dDy 
dw 



2 

a s 



-2 



aa% 



v r (Q. K r - 2v r ) 



+ 2 



aa s 

2(t-j)(Cl K r-2v r ) 
rv r a s 



y-i 

y + 1 



21(1 -j) 



dr 



v r \y + 1 



da s 
~dF 



a s r 3 



y + 1 



(23) 



where d{D.Kr)/dr = -^(r + 2)/[2(r - 2)]. The values of vy, t and a s at the viscous radius are 
determined to satisfy D v — N v — 0. 

It is noted that by using the method described above we do not use the outer boundary con- 
ditions for r > r v and the inner boundary conditions for r < r s . If we want to the solutions for 
specified outer or inner boundary conditions, we should choose the values of r s and j to satisfy the 
specified outer or inner boundary conditions. The correspondences between the values r$ and j and 
the outer solutions are investigated by Lu, Gu and Yuan (1999) for the special solutions satisfying 
the condition I = j at r = r s , and they provide solutions for the ADAF-thick disk solution, the 
AD AF- thin disk solution, and the alpha-type solution connecting the inner of outer regions. The 
special condition I = j at r = r s simplify the equations described in the last section. In the case 
of the solutions including the cases for t + j at r — r s which can be calculated by the method 
described above, we expect the basic picture is same as the solutions with £ = j at r = r s in Lu, Gu 
& Yuan (1999). 

In Fig. [TJ we present the sample solutions for ADAF-thick disk (thick lines) and ADAF-thin 
disk (thin lines) by solid lines for a — 0.1 and y = 1.5; the radial velocity v r (top-left panel), the 
angular momentum (top-right panel), the sound speed a s (bottom-left panel ) and H/r (bottom- 
right panel). The filled circles and triangles denote the positions of the sonic radius, r s , and the 
viscous radius, r v for all panels. We give the solutions in the range r/, < r < 10 1 5 , and the adopted 
values of log r s and j are log r s = 0.7 and j = 2.0. Here, we adopt the Runge-Kutta-Fehlberg 
algorithm for numerical calculations of the coupled differential equations of v r , ( and a s . The line 
of log |v r ] crosses the line of log[2y/(y + l)] l ^ 2 a s at the sonic radius, r s , and the line of logc v at the 
viscous radius, r v , respectively. In the outer region r > r v , the angular momentum of the solution 
achieves the Keplerian value. According to the all parameter search of the ADAF solutions by 
Lu, Gu & Yuan (1999), the ADAF solutions are classified into the ADAF-thick disk solution, the 
ADAF-thin disk solution and the <?-type solution connecting the inner region or the outer region. 
While the parameter spaces of r s and j for the ADAF-thick solutions and the cr-type solutions are 
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Fig. 1. The sample solutions for ADAF-thick disk (thick lines) and ADAF-thin disk (thin lines) by 
solid lines for a = 0.1 and y = 1.5; the radial velocity v r (top-left panel), the angular momentum 
(top-right panel), the sound speed a s (bottom panel ) and H/r (bottom-right panel). The filled 
circles and triangles denote the positions of the sonic radius, r s , and the viscous radius, r v . In the 
panel of log C, we also plot the line of the Keplerian angular momentum Ck(= Qur 2 ) by a dotted 
line. 

relatively wide, which are easily confirmed numerically, the parameter spaces for the ADAF-thin 
disk solutions are quite limited. By the calculation procedures presented in this paper, we can cover 
the all parameter spaces of the transonic solutions of ADAFs. 

5. Discussion 

The formula in the present study become the formula for the acausal viscous prescription if we set 
the diffusion timescale t v to be null. So, the analytic expansion around the singular points presented 
in this paper are also used to calculate the global solution for ADAF with the acausal viscosity. In 
the case of the acausal viscosity, since the boundary condition at the viscous point can not be 
used, the boundary condition at the other radius is required, e.g. no-torque condition at the horizon. 
Although the finite diffusion velocity is physically motivated, the resultant solutions have generally 
non-zero torque at the horizon which seems to be incorrect from the point of view of the relativity. 
Also, in the causal viscous prescription, the differential equation for the angular momentum become 
more complex than that of the acausal viscosity. So, although the finite diffusion speed is more 
reasonable than the infinite diffusion speed, it is noted that the calculations of the ADAF structure 
by using the causal viscous prescription may not be clearly better than the calculations by using 
the acausal viscosity. 

In the fully general relativistic calculations performed by Gammie & Popham (1998) which use 
the causal viscosity prescription, the boundary condition at the viscous point is used when calculat- 
ing the global solution. In equation (60) of Gammie & Popham (1998), although the viscous stress 
tensor at the horizon is finite if measured in the local rest frame, i.e. f( r ) ((W + 0, the viscous stress 
tensor at the horizon become null measured in the Boyer-Lindquist frame, i.e. f. = 0, because of 
T> — at the horizon. Since the location of the horizon is not the special location from the point of 
view of the local observer moving along the fluid's motion, non-zero stress in the local rest frame 
seems to be quite natural. The zero-torque condition arises only when we use the coordinate system 
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having the coordinate singularity at the horizon such as the Boyer-Lindquist coordinate and mea- 
sure the viscous stress in such frames. Then, the problem of the non-zero torque or zero-torque at 
the horizon is the issues which require the fully relativistic treatments with the spacetime structure 
and the frame-transformation between the observers measuring the viscous stress. So, it may be 
impossible to clearly resolve the torque problem at the horizon by using the calculations with the 
pseudo-Newtonian potential like the present study. 

6. Concluding Remarks 

We present the basic equations and sample solutions for the steady-state global solutions of the 
ADAFs using the causal viscous prescription and the Paczynskii-Wiita potential by the algorithm 
of explicit numerical integrations, such as the Runge-Kutta method. The calculation procedures to 
obtain these solutions are also presented. In this calculation procedure, we first solve the physical 
values at the sonic radius where L'Hopital's rule is used. The method presented in this paper enables 
us to stably solve the global solutions of ADAFs by the Runge-Kutta method, and the all parameter 
spaces of r s and j for the transonic solutions of ADAFs can be covered in this method. If we set the 
diffusion timescale to be null, the formalism in this study includes the case of the acausal viscosity 
which is usually used in the past study. Since the calculation methods in this study use the analytic 
expansion around the singular point and the numerical integration are performed by the explicit 
integration, the numerical calculations become faster. 
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